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IS.  ROSTRACT 


A  statistical  analysis  of  Alouette  topside  ionograms  '  hows  that  the  prob¬ 
ability  of  occurrence  of  the  midlatitude  F-layer  trough  ' ^  100  percent  in  win¬ 
ter  and  at  equinox,  and  is  SO  percent  in  summer.  T>ne  trough,  which  is  gen¬ 
erally  centered  on  55°  corrected  geomagnetic  latitude  and  which  is  about  5° 
to  10°  wide,  extends  throughout  the  night  hemisphere  with  f0F2  values  less 
than  2  MHz. 

Comparison  of  the  location  of  the  trough  to  the  area  in  which  coverage  by 
an  over  the  horizon  backscatter  (OTHB)  system  is  anticipated  shows  that  the 
trough  will  adversely  affect  system  performance  in  the  northeast  direction 
most  of  the  night.  An  in-depth  case  study  of  the  effects  of  the  trough  on  HF 
propagation  was  made,  based  on  an  instantaneous  three-dimensional  descrip¬ 
tion  of  the  ionosphere.  This  description  of  the  ionosphere  was  composed  of 
cross  sections  showing  ionospheric  electron-density  structure  in  the  area 
northeast  of  the  Polar  Fox  II  site.  Three-dimensional  raytracing  was  applied 
to  this  description  of  the  ionosphere,  and  the  trough  and  the  steep  electron- 
density  gradient  at  the  poleward  edge  of  the  trough  were  found  to  adversely 
affect  propagation  of  HF  signals.  For  5  MHz,  deviations  from  the  initial 
azimuth  were  found  to  be  as  big  as  15°  '(or  400km  distance  to  the  great  circle 
of  initial  azimuth)  for  rays  landing  at  the  skip  distance.  The  deviations  de¬ 
crease  with  increase  in  frequency  to  8  or  9  MHz,  but  here  propagation  de¬ 
pends  on  rays  with  low  elevation  angles  (<5°). 
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Detailed  Specification  of  the  Arctic  Ionosphere  and 
an  Application  to  Three-Dimensional  Raytracing 


1.  INTRODUCTION 

Some  OTH  backscatter  systems  in  the  planning  stage  or  in  an  experimental 

test  phase  deal  with  a  propagation  environment  quite  different  from  that  found  in 

moderate  zones-namely,  the  arctic  ionosphere  with  all  its  irregular  structure  and 

rapidly  varying  behavior.  Airborne  Investigation  by  AFCRL  has  shown  that  the 

arctic  ionosphere,  in  spite  of  its  complexity,  exhibits  a  large-scale  organization 

of  a  variety  of  ionospheric  features  that  are  either  collocated  with  or  parallel  to 

the  auroral  oval.  The  auroral-oval  concept  provides  the  framework  for  ordering 

arctic  ionospheric  phenomena.  In  fact,  the  organization  of  the  arctic  ionosphere 

is  now  well  enough  understood  that  an  analog  model  of  the  arctic  ionosphere  has 
1  2 

been  developed.  ' 

3 

Figure  1  shows  an  electron-density  cross  section,  constructed  by  Gassmann, 
depicting  the  present,  generalized  concept  of  the  structure  of  the  arctic  ionosphere 


(Received  for  publication  26  November  1973) 

1.  Gassmann,  G.  J.  (1972)  An  Ionospheric  Model  for  the  Arctic.  Environ¬ 
mental  Research  Paper,  No.  332,  APcPL-Vo-O^ez. 

2.  Gassmann,  G.  J.  (1973)  Analog  Model  1972  of  the  Arctic  Ionosphere,  pre¬ 
print,  to  be  published  as  an  AFCHL  Report  (February  1973). 

3.  Gassmann,  G.  J.  (1972)  High-latitude  ionosphere,  McGraw-Hill  Yearbook 
of  Science  and  Technology,  McGraw-Hill  Book  Co. ,  New  York,  pp,  Z57-Z60. 
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Figure  1.  Noon- To- Midnight  Cross  Section  of  the  Arctic  Ionosphere 


at  solar  maximum.  The  cross  section  applies  to  the  noon-midnight  meridian  plane 
along  approximately  80°W  and  100'E.  Along  this  meridian  plane,  in  midwinter, 
the  noon  sector  of  the  auroral  oval  is  in  darkness  due  to  the  tilt  of  the  magnetic 
axis  of  the  earth  with  respect  to  the  rotational  axis  of  the  earth.  This  cross  sec¬ 
tion  shows  circumpolar  continuity  of  certain  ionospheric  features,  for  example, 

4  5 

the  auroral  E  layer  and  the  midlatitude  F-layer  trough.  These  features  are  more 
pronounced  in  darkness  than  in  sunlight  when  they  are  obscured  by  solar-produced 
ionization. 

In  this  cross  section,  starting  from  the  right,  the  end  of  the  midlatitude  night¬ 
time  F  layer  is  seen,  sharply  bounded  in  the  north  by  an  approximately  10° -wide 
trough  in  which  fQF2  values  are  typically  <  2  MHz.  The  trough  is  terminated  in 


4.  Wagner,  R.  A.,  and  Pike,  C.  P.  (1972)  A  discussion  of  arctic  ionograms, 
AGARD  Conf.  Proc.  97,  Frihagen,  edj  pp.  4-1  to  4-20. 

5.  Muldrew,  D.B.  (1965)  F-layer  ionization  troughs  deduced  from  Alouette 
data,  J.Geophys.  Res.  70:2635-2650. 
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the  north  by  the  plasma-ring  F-layer  ionosphere  and  by  the  auroral  E  layer. 

These  features  are  produced  by  a  spatially  continuous  precipitation  of  0. 1  -  to  3  -keV 

electrons.  The  typical  critical  frequency  of  the  plasma  ring  in  the  vicinity  of  the 

poleward  trough  wall  is  4  MHz.  For  the  auroral  E  layer  extending  north  from  the 

poleward  trough  wall,  3  MHz  fQE  is  typical.  Several  degrees  north  of  the  poleward 

trough  wall,  F-region  densities  become  somewhat  enhanced,  with  fQF2  values  of 

6  to  7  MHz.  Large-scale  irregularity  structure  is  present,  which  gives  this  region 

its  name:  F-layer  irregularity  zone  (FLIZ). '  North  of  the  FLIZ,  the  polar-cap 

fQF2  values  are  approximately  2  MHz.  The  polar  cap  exhibits  very  strong  but 

rather  sporadic  E-  and  F-region  irregularities,  which  are  probably  associated 

with  the  occurrence  of  polar- cap  aurora.  At  times  the  polar  cap  has  low  electron- 

density  levels;  hence,  a  polar  cavity.  In  the  noon  sector,  equatorward  of  the  polar 

g 

cap,  there  is  a  north-to-south  sequence  of  F-,  E-,  and  D-region  ionization  bands 
produced  by  particle  precipitation.  If  the  noon  sector  is  In  darkness,  as  shown 
here,  then  the  small  depletion  of  electron  density  in  this  area  indicates  the  location 
of  the  trough.  Finally,  the  picture  is  terminated  on  the  left  by  the  northern  edge  of 
the  noontime  midlatitude  F  region. 


2.  STATISTICAL  STUDY 

Airborne  observations  have  determined  the  general  morphology  of  the  F-region 
4  9  10 

features.  '  '  To  study  in  detail  the  dependence  of  these  features  on  Kp,  UT, 
season,  solar  cycle,  and  local  time,  the  large  data  base  resulting  from  the  Alou- 
ette  topside  ionogram  analysis  was  statistically  investigated.  The  topside  iono- 
grams  were  routinely  scaled  for  fQF2  values  by  the  Communications  Research 
Centre.  Ottawa,  Canada.  The  fQF2  data  from  1962  to  1968  are  contained  on  mag¬ 
netic  tapes,  and  CRC  kindly  provided  copies  of  these  tapes  to  AFCRL  for  special 
processing.  This  processing  basically  consisted  of  transforming  the  geographic 


6.  Andrews,  M.  K. ,  and  Thomas,  J.  O.  (1969)  Electron  density  distribution 
above  the  winter  pole.  Nature  221:223. 

7.  Pike,  C.  P.  (1971)  A  latitudinal  survey  of  the  daytime  polar  F-layers, 

J.  Geophys.  Res.  76(No.  31):7745-7553. 

8.  Pike,  C.  P.  (1972)  Equatorward  shift  of  the  polar  F-layer  irregularity 
zone  as  a  function  of  the  Kp  index,  J.  Geophys,  Res.  77(No.  34):691 1-6915. 

9.  Whalen,  J.  A.,  Buchau,  J. ,  and  Wagner,  R.  A.  (1971)  Airborne  iono¬ 
spheric  and  optical  measurements  of  noontime  aurora,  J.  Atmos.  Terr.  Phys.  3.3 
(No.4):661-678. 

10.  Buchau,  J. ,  Gassman,  G.  J. ,  Pike,  C.  P. ,  Wagner,  R.A. ,  and  Whalen,  J.A. 
(1972)  Precipitation  patterns  in  the  arctic  ionosphere  determined  from  airborne 
observations,  Ann,  Geophys.  (Paris)  2Jj(No.  2):443-453. 
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latitude  and  longitude  of  the  satellite's  position  to  corrected  geomagnetic  latitude 

11  12 

(CGL)  and  corrected  geomagnetic  local  time  (CGLT),  *  enabling  a  statistical 
handling  of  the  Alouette  f0F2  values  in  corrected  geomagnetic  coordinates  rather 
than  in  geographic  coordinates.  This  approach  was  selected  because  the  high- 
latitude  ionosphere  is  geomagnetically  ordered.  The  geographic  coordinates  where 
the  Alouette  satellite  recorded  an  ionogram  were  transformed  to  CGL  and  CGLT. 
The  fQF2  data  were  then  divided  into  categories  that  should  reflect  the  various 
dependencies.  Alouette  data  from  both  northern  and  southern  hemispheres  were 
included  in  this  study.  Median  fQF2  values  were  taken  over  a  5°  range  of  CGL 
and  a  1-hr  range  of  CGLT.  Isocontours  of  f0F2  median  values  were  then  drawn, 
and  Figure  2  is  a  representative  example  of  these  isocontour  maps.  The  basic 
map  grid  is  CGL  and  CGLT.  The  contour  lines  are  in  1-MHz  increments  except 
for  the  1.  5-MHz  contour.  The  hatched  area  of  median  f0F2  <  2  MHz  is  the  statis¬ 
tical  location  of  the  trough  with  a  core  of  f0F2  <  1.5  MHz.  The  center  of  the 
trough  is  found  between  55"  and  60*  CGL,  moving  slightly  south  to  55*  for 
Kp  >  1+.  The  trough  width  is  generally  greater  than  10*.  The  polar  cavity  (low 
f0F2  values)  is  seen  in  the  vicinity  of  the  magnetic  pole.  The  two  ionization 
tongues,  which  reach  from  evening  and  from  morning  toward  the  midnight  sector, 
are  the  statistical  average  of  the  FLIZ.  Due  to  the  heavy  spread-F  conditions  in 
the  ''LIZ,  the  actual  maximum  plasma  frequencies  of  large-scale  irregularities 
are  not  scaled,  and  the  1-MHz  enhancement  at  the  trough  wall,  noted  here,  shows 
only  that  the  fQF2  is  enhanced  poleward  of  the  trough.  Further  evidence  for  the 
existence  of  a  circumpolar  FLIZ  is  the  closed  3-MHz  contour,  located  at  75*  CGL 
in  the  0900  CGLT  sector,  and  the  3-  and  4-MHz  contours  extending  toward  the 
1800  CGLT  sector. 

Past  studies  of  the  effects  of  the  arctic  ionosphere  on  OTH  and  other  radar 

13 

systems,  conducted  by  AFCRL,  have  shown  that  strong  horizontal  electron- 
density  gradients,  found  at  the  poleward  trough  wall,  and  low  electron  densities, 
jund  in  the  trough,  strongly  affect  design  parameters  and  performance  of  any 
planned  system.  Looking  again  at  this  general  problem,  the  situation  of  an  OTHB 
system  assumed  to  be  located  at  the  POLAR  FOX  II  site  was  considered. 


11.  Hakura,  Y.  (1965)  Tables  and  maps  of  geomagnetic  coordinates  cor¬ 

rected  by  the  higher  order  spherical  harmonic  terms,  Rep.  Ionosph.  Space  Res. 
(Japan)  19:121-157.  - 

12.  Whalen,  J.  A.  (1970)  Auroral  Oval  Plotter  and  Nomograph  for  Determin¬ 
ing  Corrected  Geomagnetic  Local  Time,  Latitude,  and  Longitude  in  the  Northern' 
Hemisphere,  AFdlfrL-70-0422,  Environmental  Research  Paper  No.  327,  Bedford. 
Mass. 

13.  Evans,  D.  L.  (1971)  Comp.,  Ionospheric  and  tropospheric  limitations  to 
radar  accuracy.  Air  Force  Surveys  in  Geophysics,  No.  231,  AFCRL-71-0169. 
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Figure  2.  Contours  of  Median  fftF2  (in  MHz)  Derived  From  Alouette  1 
Topside  Soundings.  Data  base:  Winter  1963,  0000  to  1200  UT,  k„  0  to 
1+.  The  coordinate  system  is  corrected  geomagnetic  latitude  ana  cor¬ 
rected  geomagnetic  local  time 


The  radar  beam  has  been  directed  toward  the  northeast  and  has  an  azimuth  range 
of  ±  45*.  The  area  of  coverage  of  this  radar  beam  has  been  converted  into  the  CG 
coordinate  system  and  is  superimposed  on  the  fQF2  contour  map  in  Figure  2. 

This  presentation  clearly  points  out  that  the  boresight  is  aligned  with  the  center  of 
the  trough,  and  the  situation  does  not  improve  for  southerly  directions.  Trans¬ 
mission  toward  the  north  will  encounter  the  poleward  trough  wall  and  the  FLIZ, 
and  will  therefore  be  strongly  affected  by  auroral  activity. 

The  actual  situation  of  strong  north-south  electron-density  gradients  is  not 
adequately  represented  by  this  averaged  picture.  The  steepest  gradients  found  in 
the  statistical  analysis  were  about  0.5  MHz  per  100  km.  The  actual  gradients. 


14  4 

determined  from  ground-based  measurements  and  airborne  measurements  are 

15 

approximately  2  to  3  MHz  per  100  km. 

The  probability  of  occurrence  of  f0F2  values  s  2  MHz  was  then  determined 
(Figure  3).  Based  on  1963  (solar  minimum)  data  and  on  1968  (solar  maximum) 
data,  the  trough  is  present  about  90  percent  of  the  time  in  the  winter  and  at  equi¬ 
nox  with  a  100  percent  probability  of  occurrence  near  midnight.  The  probability 
of  occurrence  of  the  trough  drops  to  50  percent  in  summer.  The  local-time  extent 
of  the  trough  has  a  strong  seasonal  dependence  due  to  the  change  in  solar  illumina¬ 
tion.  The  trough  extends  for  approximately  12  hr  and  is  centered  on  01  CGLT 
in  winter,  extends  for  6  to  8  hr  at  equinox  and  is  centered  on  02  CGLT.  and 
extends  in  summer  for  less  than  6  hr  and  is  centered  on  03  CGLT.  Further 
study  of  the  variability  and  the  dependencies  of  the  trough  is  still  in  progress.  In 
this  paper,  effects  of  the  trough  on  propagation  will  be  discussed  for  a  specific 
case  study. 

The  statistical  picture  has  shown  that  any  OTHB  system  planned  to  operate 
at  these  latitudes  will  have  to  contend  with  the  trough  ionosphere.  Therefore,  an 
attempt  was  made  to  specify  a  large  portion  of  the  arctic  ionosphere  that  could 
then  be  applied  to  an  OTH  situation.  This  was  done  for  two  reasons.  First,  we 
wanted  to  show  by  a  well-documented  example  that  median  f0F2  contour  maps  need 
interpretation  because  of  the  steep  gradients  that  exist  in  reality.  Second,  we 
wanted  to  provide  the  input  data  for  a  three-dimensional  raytracing  program  and 
to  clearly  show  the  effects  of  trough  and  strong  north-south  gradients  on  OTH 
propagation. 

3.  IONOSPHERIC  SPECIFICATION 

On  9  December  1971,  a  coordinated  data -acquisition  program  involving  air¬ 
craft,  satellite,  and  ground  stations  provided  a  detailed  picture  of  the  ionosphere 
in  the  area  north  and  east  from  the  POLAR  FOX  II  radar  site.  Airborne, 
satellite-borne,  and  ground-based  ionosondes  measured  electron-density  profiles, 
and  from  these  measurements  the  trough,  the  steep  electron-density  gradient 
poleward  of  the  trough,  and  the  auroral  ionosphere  could  be  specified.  In  addition, 
unique  new  data  that  specified  auroral  conditions  over  a  vast  area  were  used. 


14.  Bowman,  G.G.  (1969)  Ionization  troughs  below  the  F2-layer  maximum. 
Planet.  Space  Sci.  ^7:777-796. 

15.  Buchau,  J.  (1972)  Instantaneous  versus  averaged  ionosphere,  from 
"Arctic  Ionosphere  Modeling-Five  Related  Papers",  Air  Force  Surveys  in  Geo- 
physics.  No.  241,  AFCRL-72-0305. 
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Figure  3.  Contours  of  Percentage  Occurrence  of  f0F2  <  2  MHz. 
Data  base  and  coordinate  system  same  as  Figure  2 


These  data  are  collected  by  n  scanning  photometer  on  board  a  satellite.  An  ex¬ 
ample  of  these  new  data  is  presented  in  Figure  4.  This  picture  covers  the  eastern 
half  of  the  United  States  and  Canada  and  a  major  part  of  the  Atlantic  Ocean,  in¬ 
cluding  Greenland  and  the  Arctic  Ocean  up  to  the  North  Pole.  Of  specific  interest 
in  the  picture  is  an  east-to-west-orlented  belt  of  continuous  aurora  extending 
from  Central  Canada  to  Greenland,  Near  the  poleward  edge  of  the  continuous 
aurora,  discrete  auroral  arcs  are  seen. 

The  band  of  continuous  aurora  results  from  excitation  of  the  arctic  atmosphere 
by  electron  precipitation.  The  ionization  that  results  from  this  precipitation  pro¬ 
duces  the  auroral-E  layer,  which  is  closely  associated  with  the  continuous  aurora 
and  the  auroral- F  layer.  The  poleward  trough  wall  is  approximately  collocated 
with  the  equatorward  edge  of  the  continuous  aurora.  In  this  picture,  moonlight 
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Figure  4.  Satellite  Picture  Produced  by  Scanning 
Photometer,  Showing  Midnight  Sector  of  Auroral 
Oval.  A  broadband  of  continuous  aurora  is  bounded 
in  the  north  and  partially  in  the  south  by  discrete 
aurora.  Geographic  features  indicate  dimensions 
of  the  coverage  area 
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reflecting  from  clouds  occurs  in  the  east.  By  relating  the  location  of  auroral 
structure  present  in  the  picture  to  electron-density  levels  and  gradients  mea¬ 
sured  by  ionosondes  and  the  ISIS  2  topside  sounder,  a  large-scale  and  accurate 
specification  of  the  ionosphere  could  be  inferred  over  the  area  of  the  picture  in 
which  no  ionosonde  data  were  recorded.  In  addition,  by  relating  the  location  of 
auroral  structures  to  the  location  of  the  radar  site,  the  exact  range  from  the  site 
to  the  location  where  auroral  effects  would  be  expected  to  affect  propagation  con¬ 
ditions  was  read  directly  from  the  picture.  Hence,  a  specification  of  the  iono¬ 
sphere  and  of  HF  propagation  conditions  of  unprecedented  dimensions  has  been 
achieved.  These  large-scale  pictures  of  auroral  conditions  should,  in  general,  be 
made  available  on  a  near-real-time  basis  to  any  system  that  may  be  affected  by 
the  ionosphere  or  by  the  aurora.  By  utilizing  these  pictures  and  routine  frequency 
predictions,  an  accurate  estimate  of  propagation  conditions  in  close  to  real  time 
could  be  reached. 

A  radar-beam  position,  directed  northeast  (45°  east  of  north)  from  the  POLAR 
FOX  II  site,  was  selected  as  the  center  line  of  a  76°  sector  in  which  a  three- 
dimensional  picture  of  the  ionosphere  was  constructed.  The  sector  is  noted  in 
Figure  5  as  the  area  between  the  dashed  lines.  The  northern  hemisphere  has 
been  drawn  in  corrected  geomagnetic  coordinates.  The  central-beam  position  is 
labeled  0°,  and  beam  positions  ±  38.0“  from  the  central-beam  position  are  also 
noted. 

Shading  corresponds  to  the  location  of  the  continuous  aurora  seen  in  Figure  4. 
In  the  azimuth  sector  -6“  to  +6“  with  respect  to  the  selected  boresight  (45°  east 
of  north),  electron-density  cross  sections  were  provided  in  2°  increments  (seven 
cross  sections).  The  remaining  sector  between  -38“  and  +38°  was  described 
every  4“  in  azimuth.  Selected  electron-density  cross  sections  representative  of 
the  total  of  25  cross  sections  produced  are  given  in  Figure  6.  The  -34°  cross 
section  (top  of  the  figure)  is  a  south-to-north  cut  through  the  polar  ionosphere 
similar  to  the  generalized  presentation  shown  in  Figure  1.  The  -4°  and  0“  cross 
sections  just  skim  the  poleward  trough  wall,  showing  sections  of  the  trough  iono¬ 
sphere  to  the  south  of  it. 

The  ionosphere  south  of  the  poleward  trough  wall  is  shown  in  the  +34°  cross 
section.  It  is  in  general  a  2-MHz  ionosphere  with  a  narrow  region  of  1.  5-MHz 
f0F2.  This  trough  ionosphere  fluctuates  in  height  between  360  and  470  km,  which 
are  the  values  taken  from  the  ISIS  2  topside  soundings  of  the  0420  UT  pass  that 
went  exactly  over  the  POLAR  FOX  II  site,  the  origin  of  these  cross  sections. 

In  these  cross  sections  the  poleward  trough  wall  is  characterized  by  an  abrupt 
increase  in  plasma  frequency  from  2  MHz  in  the  trough  to  4  MHz  in  the  auroral 
F-layer  or  plasma  ring.  The  FLIZ  is  seen  in  the  -34“  cross  section  as  a  5-MHz 
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Figure  5.  Continuous  Aurora  Shown  in  Figure  4  Transferred  into  the  Corrected 
Geomagnetic  Coordinate  System.  Note  extremely  good  representation  of  tha 
aurora  midnight  sector  in  corrected  geomagnetic  coordinates.  Black  square  in¬ 
dicates  location  of  POLAR  FOX  II.  The  sector  of  the  ionosphere  studied  lies 
between  -38°  and  +38“  azimuth  with  respect  to  the  0°  boresight  in  25  radial 


cross  sections 
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enhancement  north  of  the  aurora.  The  aurora  itself  is  represented  as  an  auroral 
Es  with  ftEga  about  3  to  4  MHz  at  1800  km  distance  from  the  POLAR  FOX  II 
site  in  the  -34*  azimuth. 

The  strong  horizontal  gradients  of  2  to  3  MHz  per  100  km  found  in  the  vicin- 

14  4 

ity  of  the  poleward  trough  wall  by  Bowman  and  Wagner  and  Pike  are  definitely 
present  in  this  case  study  of  9  December  1971.  In  order  to  study  the  effects  on 
radio  wave  propagation  of  this  specific,  but  (from  our  experience,  for  this  loca¬ 
tion  and  time)  rather  typical  ionosphere,  the  cross  sections  were  used  as  an  input 
to  a  three-dimensional  raytracing  program. 


4.  RAYTRACING 

In  a  first  raytracing  attempt  a  spherical- coordinate  system  was  positioned 
with  its  pole  at  the  radar  site.  Seven  ionospheric  cross  secii-ns  (±  6°,  ±  4s,  ±  2s, 
and  0*  with  respect  to  a  boreslght  of  45*  east  of  north)  were  used  as  the  meridi¬ 
onal  planes  of  this  coordinate  system.  Since  it  turned  out  that  the  deviation  of  the 
computed  rays  was  greater  than  the  described  ionosphere  allowed  for,  the  iono¬ 
sphere  was  described  by  a  total  of  25  cross  sections  in  an  azimuth  range  *  38* 
from  boreslght.  This  ionosphere  was  transformed  into  the  geomagnetic  dipole 
system  for  a  second  raytracing  attempt.  Both  chosen  coordinate  systems  are 
spherical,  so  that  all  parts  except  the  electron-density  subroutine  of  the  highly 
regarded  Jones  raytracing  program18  can  be  used  for  our  purpose.  The  dipolar 
coordinates  are  not  as  deuirable  as  the  corrected  geomagnetic  coordinates,  but 
are  a  considerable  imprc  vement  over  the  geographic  coordinates  for  simulation 
of  electron  distributions  in  the  polar  ionosphere.  In  particular  the  boundary  of  the 
trough  is  nearly  constant  in  corrected  geomagnetic  latitude,  and  extends  over  an 
Interval  of  several  degrees  in  dipolar  latitude  but  tens  of  degrees  in  geographic 
latitude.  The  selection  of  a  coordinate  system  that  better  orders  the  ionospheric 
features  simplifies  and  Improves  the  interpolation  between  two  meridional  planes. 

AFCRL  has  developed  an  electron-density  subprogram  and  Incorporated  it 
into  the  Jones  raytracing  program.  Input  data  to  this  subprogram  consist  of 
plasma-frequency  values  specified  at  uniform  grid-size  intervals  of  r,  9,  and  <t> ; 
in  the  present  etudy.  Intervals  of  4  km  in  height,  0.5*  dipolar  latitude,  and  5*  in 
longitude  were  used.  The  output  of  this  subprogram,  feeding  into  the  Jones  ray¬ 
tracing  program,  provides  values  of  the  plasma  frequency  and  its  spatial  deriva¬ 
tives  at  fine-step  Intervals  (determined  by  Integration  step  size)  along  the  raypath. 


16.  Jones,  R.  M.  (1968)  A  three-dimensional  raytracing  computer  program. 
Radio  Science  3(No.  l):93-94. 
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The  input  to  the  electron-density  subprogram  comes  from  an  external  pre¬ 
processing  program.  The  input  to  the  preprocessing  program  consists  of  the  loci 
of  plasma-frequency  contours  in  each  of  a  number  of  vertical  meridional  planes 
of  the  coordinate  system  used.  Meridional  planes,  and  data  points  in  each  plasma- 
frequency  contour,  are  chosen  at  all  locations  where  significant  changes  occur  in 
the  principal  features  (such  as  the  trough  and  plasma  ring)  of  the  electron  distri¬ 
bution. 

Within  the  preprocessing  program,  a  particular  mode  of  interpolation  intro- 
17 

duced  at  AFCRL  is  used  in  order  to  ensure  that  a  principal  feature  that  appears 
in  each  of  two  chosen  meridional  planes  would  be  projected  as  one  and  only  one 
principal  feature  at  all  intermediate  planes  between  the  two  chosen  ones.  If  this 
stipulation  were  not  fulfilled,  incorrect  horizontal  gradients  of  the  electron-density 
(hence  incorrect  off-great- circle  deviation  of  raypaths)  would  be  obtained  at  some 
portions  of  the  raypaths.  Within  the  electron-density  subprogram  of  the  raytracing 
program,  ordinary  interpolation  is  used.  Hence,  at  the  level  of  fine  coordinate 
increments  (4  km  in  height,  0.5°  in  0,  and  5*  in  0  for  the  present  study),  residual 
and  extraneous  errors  in  horizontal  gradients  of  electron-density  still  occur;  how¬ 
ever,  these  residual  errors  would  cause  only  random  fluctuations  in  the  off-great- 
circle  deviations  of  raypaths. 

The  principal  drawback  of  the  methoJ  of  simulating  electron-density  distribu¬ 
tions  described  above  is  the  large  size  of  the  central  memory  required  in  the  cen¬ 
tral  processor  of  the  computer  system.  A  central-memory  size  of  260, 000  (in 
the  octal-number  system)  was  used,  approximately  200,  000  for  numerical  data, 
and  50,  000  for  the  Jones  raytracing  program. 

There  are  two  chief  advantages  of  the  method:  raytracing  speed,  and  a  do-it- 
yourself  capability  for  the  raytracing  user.  It  takes  approximately  0.7  sec  of 
computing  time  with  the  CDC-6600  computer  to  compute  a  one-hop,  3000-km 
raypath.  The  user  can  easily  visualize  the  simulated  electron-density  distribution 
that  is  actually  introduced  into  the  raytracing  computations,  and  retains  assurance 
that  the  electron-density  would  be  well  approximated  at  the  data-grid  points  and 
also  that  there  would  be  no  wild  or  unintended  fluctuations  of  electron-density 
between  grid  points.  The  so-called  problem  of  control  or  instability,  troublesome 
in  analytical  representations  of  electron-density  involving  summation  of  a  large 
number  of  terms  (each  a  product  of  functions  of  r,  0,  and  0),  is  bypassed  in  the 
method  used  here. 


17.  Wong,  M.  S.  (1973)  AFCRL  Report  to  be  published. 
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5.  RAYTRACING  RESULTS 


Figure  7  shows  as  an  example  the  results  of  the  raytracing  for  6  MHz,  start¬ 
ing  with  an  initial  azimuth  of  -1*  relative  to  the  cho3en  boresight  of  45°  east  of 
north.  The  upper  part  shows  a  projection  of  the  raypaths  into  the  meridional 
plane  of  the  initial  azimuth  vhile  the  lower  part  shows  the  deviations  from  the 
initial  azimuth  in  a  projection  into  the  horizontal  plane,  given  here  in  range  and 
azimuth.  A  group  of  curves  allows  one  to  judge  the  distance  of  the  landing  points 
of  a  ray  to  the  great  circle  of  the  initial  azimuth.  The  following  figures  will  show 
only  the  heavy  line  connecting  the  landing  points  in  Figure  7  as  the  most  rele¬ 
vant  result  of  the  raytracing.  Raytracing  was  done  for  5,  6,  7,  8,  and  9  MHz 
for  initial  azimuths  of  +5  to  -9“  (in  1®  increments)  with  respect  to  the  0®  refer¬ 
ence  (45®  east  of  north). 
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Figure  7.  Example  of  Three-Dimensional  Raytracing.  Top  shows  projec 
tion  of  raypaths  into  the  vertical  plane  of  Initial  azimuth.  Bottom  shows 
projection  into  the  horizontal  plane 
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Figure  8  gives  the  results  for  ti  MHz.  The  numbers  at  the  right  ends  of  the 
loci  of  the  landing  points  indicate  the  initial  azimuth  with  respect  to  45°  east  of 
north.  The  deviations  are,  as  expected,  mostly  toward  the  south  of  the  initial 
azimuth  and  are  largest  for  the  rays  launched  most  northerly  and  are  decreasing 
for  rays  intersecting  less  and  less  of  the  enhanced  F  region  north  of  the  trough 
wall.  A  test  ray  was  launched  at  +13°  azimuth  so  as  to  encounter  only  the  rather 
horizontally  strati. led  2  MHz  trough  ionosphere.  As  expected,  there  were  no 
deviations,  and  the  range  covered  extended  from  2200  to  3900  km,  as  indicated 
by  the  heavy  line  on  the  range  axis.  Lauding  poi'its  for  rays  -9#  to  -7°  fall  into 
two  branches,  the  big  looped  curves  resulting  from  rays  leaving  under  low  eleva- 
ticn  angles  (0°  to  10°),  and  the  long  curves  resulting  from  rays  with  elevation  of 
12®  to  25®.  Rays  between  10®  and  12®  elevation  are  M-mode  propagation. 


Figure  8.  Raytracing  Results  for  5  MHz.  Shown  are  the  loci  of  the 
landing  points  as  a  function  of  the  initial  azimuths.  Initial  azimuths  with 
respect  to  a  central  boresight  of  45®  east  of  north  are  Indicated  by  the 
numbers  to  the  right  of  the  curves 
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skipping  the  earth  due  to  ionospheric  tilts  at  the  reflection  point.  The  skip  dis¬ 
tance  can  clearly  be  seen;  the  hooks  at  the  upper  ends  of  the  landing-point  curves 
result  from  *he  unfolding  of  the  Pedersen  ray  in  this  diverging  ionosphere.  For 
rays  -6°  t  >  -1°,  all  rays  with  elevations  smaller  than  about  10°  are  M  modes, 
the  landing  points  for  the  higher  rays  forming  a  continuous  curve.  For  azimuths 
08  to  5°,  low  rays  (0°  to  5°  elevation)  reach  the  ground  again,  high  rays  (14°  to 
20°  <  _evation)  giving  somewhat  poor  close-in  coverage. 

At  6  MHz  (Figure  9),  the  situation  is  the  same  in  general.  The  skip  distance 
has  increased,  azimuths  -9°,  -8°.  and  -7“  result  again  in  two  branches,  one 
branch  from  low  rays  and  a  second  branch  from  higher  rays,  the  break  being  due 
to  M-mode  rays.  Rays  of  -6°  to  -l8  azimuth  land  only  for  elevations  greater 
than  about  58.  Rays  of  0°  to  5°  azimuth  have  again  a  branch  due  to  low  elevations 
(about  08  to  8°),  which  at  times  are  connected  to  the  branch  of  higher  elevations. 
The  test  ray  propagated  at  13“  into  the  horizontally  stratified  trough  shows  no 
deviation  and  covers  tne  range  from  3100  to  4300  km. 
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Figure  9.  Raytracing  Results  for  6  MHz 
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Figure  10  shows  the  ground  range  and  deviation  for  7  MHz.  Since  this  fre¬ 
quency  is  close  to  the  MUF,  high  rays  produce  only  broken-up  coverage;  for  rays 
of  -9°  to  -6“  azimuth,  the  coverage  is  provided  by  low  rays  (0°  to  5°  elevation), 
with  no  high  ray  reaching  the  ground.  The  2-MHz  trough  can  no  longer  support 
propagation  at  this  frequency. 
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The  8-MHz  results  (Figure  11)  show  that  only  rays  of  -9#  to  -6°  azimuth  are 
reflected,  the  deviation  decreasing  rapidly.  The  coverage  is  provided  by  elevation 
angles  0s  to  5°.  Finally,  at  9  MHz  (Figure  12),  skip  distance  is  almost  3000  km; 
only  rays  of  -9°  to  -5°  azimuth  are  reflected,  relying  on  elevation  angles  <5°. 
For  10  MHz  no  reflection  is  achieved  for  the  assumed  initial  azimuth  of  -9®  to 
+  13°  with  respect  to  45®  east  of  north. 
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GROUNDRANGE  AND  DEVIATION  FROM  GREAT-CIRCLE  PATH 
1000 _ 2000 _ 30£0 _ 4000  km 


6.  CONCLUSIONS 

The  Ne-gradient  at  the  poleward  trough  wall  haB  been  shown  to  severely 
affect  the  propagation  of  radio  waves.  Deviations  from  the  original  azimuth  are 
shown  to  be  as  big  as  25°  or  1200  km  off-great-circle  for  5  and  6  MHz.  For  the 
rays  landing  at  the  skip  distance,  deviations  are  as  big  as  15s  or  400  km  off- 
great- circle.  For  8  and  9  MHz  the  deviations  are  smaller,  but  only  rays  with 
low  elevation  angles  propagate  to  the  ground.  It  has  to  be  pointed  out,  however, 
that  the  azimuths  selected  present  the  worst  case  as  far  as  the  deviation  is  con¬ 
cerned.  Since  the  trough  is  a  routinely  occurring  and  variable  feature  of  the 
southern  edge  of  the  polar  ionosphere,  an  operational  system  requires  instant  cali¬ 
bration.  This  in  turn  requires  either  a  large  number  of  transponders  or  a  real¬ 
time  determination  of  the  location  of  the  trough  together  with  a  library  of  precal¬ 
culated  results  from  raypaths  such  as  those  shown  here.  A  model  of  the  arctic 
ionosphere  and  methods  to  determine  the  specifics  in  real  time  are  being  worked 
on  at  AFCRL. 
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DEVIATION  FROM  INITIAL  AZIMUTH 


The  trough  ionosphere  supports  propagation  at  5  and  6  MHz  with  skip  dis¬ 
tances  of  2000  to  3000  km,  clearly  demanding  a  high  emphasis  of  the  low-frequency 
capability  of  any  system  planned  to  operate  in  this  region. 

The  deviations  described  result  in  double  coverage  of  certain  areas  south  of 
the  trough  wall  and  in  regions  not  covered  at  all  in  the  direction  of  the  modeled 
beams. 

It  is  planned  to  complete  the  computations  for  all  azimuths,  to  determine  total 
coverage  possible  for  'he  available  test  case.  POLAR  FOX  II  data  are  available 
for  9  December  197 1,  and  a  comparison  of  the  computer  simulation  and  the  actual 
data  should  show  how  well  the  modeling  depicted  the  actual  situation. 

The  computer  program  used  allows  the  input  of  an  arbitrary  ionosphere  rather 
than  a  mathematical  model  with  all  its  shortcomings,  and  after  the  presently  very 
high  requirements  on  core  memory  have  been  reduced  by  proper  use  of  disc  stor¬ 
age  it  should  be  useful  for  producing  a  comprehensive  library  of  results,  or  may 
even  be  used  for  real-time  applications. 
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Finally,  it  should  be  pointed  out  that  the  demonstrated  feasibility  of  calibrating 
an  OTH  radar  field  makes  practical  sense  because  a  reasonably  accurate  descrip¬ 
tion  of  the  ionosphere  now  appears  possible.  Of  particular  aid  are  satellite  pic¬ 
tures  showing  the  large-scale  extent  of  the  aurora.  These  pictures  are  available 
from  the  AWS,  and  their  use  directly  at  the  site  or  indirectly  after  AWS  interpre¬ 
tation  should  be  contemplated  for  any  operational  system.  The  interpretation  of 
these  pictures  in  terms  of  the  corresponding  ionospheric  features  has  been 
achieved  at  AFCRL,  and  further  refinement  of  the  interpretation  is  underway. 
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